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Abstract — This paper presents an analytical model and a 
geometric numerical integrator for a rigid body connected to an 
elastic string, acting under a gravitational potential. Since the 
point where the string is attached to the rigid body is displaced 
from the center of mass of the rigid body, there exist nonlinear 
coupling effects between the string deformation and the rigid 
body dynamics. A geometric numerical integrator, refereed to as 
a Lie group variational integrator, is developed to numerically 
preserve the Hamiltonian structure of the presented model 
and its Lie group configuration manifold. These properties are 
illustrated by a numerical simulation. 

I. Introduction 

The dynamics of a body connected to a string appear in 
several engineering problems such as cable cranes, towed 
underwater vehicles, and tethered spacecraft. It has been 
shown that gravitational forces acting along a string can alter 
the tension of the string, while significantly disturbing the 
dynamics of a body connected to the string [1]. Therefore, it 
is important to model the string dynamics accurately as well 
as the dynamics of the body even if the tension of the string 
is low. 

Several dynamic and numerical models have been devel- 
oped. Lumped mass models, where the string is spatially 
discretized into connected point masses, were developed 
in [2], [3], [4]. Finite difference methods in both the spatial 
domain and the time domain were applied in [5], [6]. Finite 
element discretizations of the weak form of the equations of 
motion were applied in [6], [7]. String deployment models 
were developed in [8], [9]. 

The goal of this paper is to develop an analytical model 
and a numerical simulation tool for a rigid body connected to 
a string acting under a gravitational potential. This dynamic 
model is referred to as a 3D elastic string pendulum: it is a 
generalization of a 3D pendulum model introduced in [10] to 
include the effects of string deformations; it is an extension 
of a string pendulum model with a point mass bob [6]. 

We assume that the point where the string is attached 
to the rigid body is displaced from the center of mass of 
the rigid body so that there exist nonlinear coupling effects 
between the string deformation dynamics and the rigid 
body dynamics. This provides a more realistic and accurate 
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dynamic model. We show that the governing equations of 
motion can be developed according to Hamilton's variational 
principle. 

The second part of this paper deals with a geometric 
numerical integrator for the 3D elastic string pendulum. Ge- 
ometric numerical integration is concerned with developing 
numerical integrators that preserve geometric features of a 
system, such as invariants, symmetry, and reversibility [11]. 
For a numerical simulation of Hamiltonian systems evolving 
on a Lie group to exhibit good long-time energy behavior, 
it is critical to preserve both the symplectic property of 
Hamiltonian flows and the Lie group structure [12]. A 
geometric numerical integrator, referred to as a Lie group 
variational integrator, has been developed for a Hamiltonian 
system on an arbitrary Lie group in [13]. 

A 3D elastic string pendulum is a Hamiltonian system, and 
its configuration manifold is expressed as the product of the 
special Euclidean group S0(3) and the space of connected 
curve segments on R'^. This paper develops a Lie group 
variational integrator for a 3D elastic string pendulum based 
on the results presented in [13]. The proposed geometric 
numerical integrator preserves symplecticity and momentum 
maps, and exhibits desirable energy conservation properties. 
It also respects the Lie group structure of the configuration 
manifold, and avoids the singularities and computational 
complexities associated with the use of local coordinates. 

In summary, this paper develops an analytical model and 
a geometric numerical integrator for a 3D elastic string pen- 
dulum. These provide a mathematical model and a reliable 
numerical simulation tool that characterizes the nonlinear 
coupling between the string dynamics and the rigid body 
dynamics accurately. This can be naturally extended to con- 
trolled dynamics, and serve as the basis for optimal control 
algorithms as in [14]. 

This paper is organized as follows. A 3D elastic string 
pendulum is described in Section An analytical model and 
a Lie group variational integrator are developed in Sectionlllll 
and in Section IIVI respectively, followed by a numerical 
example in Section FVl 

II. 3D Elastic String Pendulum 

Consider a rigid body that is attached to an elastic string. 
The other end of the string is fixed to a pivot point. We 
assume that the rigid body can freely translate and rotate in 
a three dimensional space, and the string is extensible and 
flexible. The bending stiffness of the string is not considered 
as the diameter of the string is assumed to be negligible 
compared to its length. The point where the string is attached 
to the rigid body is displaced from the center of mass of the 
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Fig. 1. 3D Elastic String Pendulum 



rigid body so that the dynamics of the rigid body is coupled 
to the string deformations and displacements. 

This model is a generalization of the 3D pendulum and 
the string pendulum introduced in [10] and [6], respectively, 
and it is referred to as a 3D elastic string pendulum. This is 
illustrated in Fig. [T] 

We choose a global reference frame and a body-fixed 
frame. The origin of the body-fixed frame is located at the 
end of the string where the string is attached to the rigid body. 
Since the string is extensible, we need to distinguish between 
the arc length for the stretched deformed configuration and 
the arc length for the unstretched reference configuration. 
Define 

/ G M Total length of the unstretched string 

s G [0, Z] Length of the string from the pivot to a ma- 

terial point P for the unstretched reference 
configuration 

s(s, i) G M+ Length of the string from the pivot to a 
material point P for the stretched deformed 
configuration 

r(s, t) G IR-^ Vector from the pivot to a material point P 
in the global reference frame 

R G S0(3) Rotation matrix from the body-fixed frame 
to the reference frame 

f2 G M"^ Angular velocity of the rigid body repre- 

sented in the body-fixed frame 

Pc G M'^ Vector from the origin of the body fixed 

frame to the mass center of the rigid body 
represented in the body fixed frame 

JI G IR+ Mass density of the string per unit un- 

stretched length 

M G M+ Mass of the rigid body 

J G R^^'^ Inertia matrix of the rigid body represented 
in the body fixed frame 

A configuration of this system can be described by the 
locations of all the material points of the string, r(s, t) for 
s G [0, 1], and the attitude of the rigid body R{t) with respect 
to the reference frame. So, the configuration manifold is 



G = C°°([0,Z],]R3) X S0(3), where C°°([0, Z], R^) denotes 
the space of smooth connected curve segments on M.^ and 
S0(3) = {i? G M^^^ I ^ J ^ ^ 

III. Continuous-time Analytical Model 

In this section, we develop continuous-time equations of 
motion for a 3D elastic string pendulum. The equations 
for a string pendulum connected to a point mass has been 
developed in [6]. Here we focus on generalizing them for 
a rigid body. The attitude kinematics equation of the rigid 
body is given by 

R = rCi, (1) 

where the hat map " : R"^ ^ so(3) is defined by the condition 
that xy = X xy for any x,y G M.^. For notational simplicity, 
we do not express the time dependency of variables explic- 
itly, i.e. r(s) = r(s, t). 

A. Lagrangian 

Kinetic energy: The total kinetic energy is composed of 
the kinetic energy of the string Tgt,. and the kinetic energy 
of the rigid body Trb- Let r(s, be the partial derivative 
of r(s, t) with respect to t. This represents the velocity of a 
material point on the string. Then, the kinetic energy of the 
string is given by 

fi 1 

Tstr = 2^ll^(s)f ds. (2) 

Let p G R"^ be the vector from the mass center of the rigid 
body to a mass element of the rigid body represented in the 
body fixed frame. The location of the mass element is given 
by r{l)+ R{pc + p) in the global reference frame. Therefore, 
the kinetic energy of the rigid body can be written as 

Trh= I \ \Hl) + R^{Pc + p)fdM{p) 

= \^Mr{l) ■ f{l) + ^n- jn + Mf{l) ■ Rnpc, (3) 

where B denotes the region enclosed by the rigid body 
surface, and we use the following properties: p dM = 0; 
xy = -yx- J = - /^((p + p,)A)2 dM. 

Potential Energy: The strain of the string at a material 
point located at r{s) is given by 

e ~ lim = s (,s) — 1, 

As^O As \ J ^ 

where ( )' denote the partial derivative with respect to s. The 
tangent vector at the material point is given by 

dr{s) dr{s) ds t'{s) 
^* ds &s ds{s) s'(s) 

Since this tangent vector has the unit length, we have s'(s) = 
II r'(s) II . Therefore, the strain is given by e = ||'''(s)|| — 1- 
The potential energy of the string is composed of the elastic 
potential and the gravitational potential; 

1 

Vstr = -EA{\\r'{s)\\ - lf-Jlgr{s) ■ eg ds, (4) 



where E and A denote the Young's modulus and the sectional 
area of the string, respectively, and the unit vector 63 
represents the gravity direction. 

Since the location of the center of mass of the rigid body 
is r{l) + Rpc in the global reference frame, the gravitational 
potential energy of the rigid body is 



Vrb = -M,9(r(0 + i?Pc) -63. 



(5) 



From (O-©, the Lagrangian of the 3D elastic string 
pendulum is given by 



T, 



rb 



Vr 



rb- 



(6) 



L — Tstr ^ Vstr 

B. Elder- Lag range Equations 

Let the action integral be © = J^^^ Ldt. It is composed 
of two parts, ©sir and ©rf,, contributed by the string and 
by the rigid body, respectively. According to the Hamilton's 
principle, the variation of the action integral is equal to 
zero for fixed boundary conditions, which yields the Euler- 
Lagrange equations of the 3D elastic string pendulum. 

By repeatedly applying integration by parts, the variation 
of ©str can be written as 



5^s 



-EA 



V{i)\ 



1 



63 



EA 



\\r'{l)\ 

\r'm - 

\\r'[s)\\ 



r'{l)-5r{l) 



I , 



1 



r'(s) ■5r{s)dsdt. (7) 



(See [6] for details.) 

Next, we found the variation of ©^b. It can be written as 



Mr {I) + MRClpc 



6r{l) + Mges ■ 

+ [.m + MpcR^r{i)] ■ sn 

+ Mf{l) ■ SRnpc + Mge3 ■ SRp^ dt. 
The variation of a rotation matrix can be written as 



5r{l) 



(8) 



5R 



d 



R' 



!=0 



d_ 

It 



R cxp ef) ~ Rfj 



c=0 



for ?7 G [15]. The corresponding variation of the angular 
velocity is obtained from the kinematics equation ([T]l: 



6n = ^ 

de 



(R'fR' = {fi + Qxr]y 



e=0 



Substituting these into ([Sj and applying the integration by 
parts, we obtain 



variation according to Hamilton's principle. This yields the 
following Euler-Lagrange equations: 



p,r{s, t) — p,ge3 — EA 



M[r{l,t)-Rpc(l 



EA' 



d_ (\\r'{s,i 
ds 



1 



r'{s,t) 



rCi^Pc 



563 

1 



lk'(M)ll 



0, 

(10) 
(11) 



jn + O JO + mpcR r{l, t) 



mgpcR^es 



0. (12) 



Conserved quantities: The total energy, given hy E = 
str + Vstr + Trb + Vrb, is preserved. As the Lagrangian is 
invariant under the rotation about the gravity direction, the 
total angular momentum about the gravity direction is con- 
served. It is given by tts = {/p jLf{s)r{s) ds + Mf{l){f{l) + 
RClp^) - Mi-{l)Rp^ + R.m} ■ 63. 

IV. Lie Group Variational Integrator 

The continuous-time Euler-Lagrange equations developed 
in the previous section provide an analytical model for a 
3D elastic string pendulum. However, the popular finite dif- 
ference approximations or finite element approximations of 
those equations using a general purpose numerical integrator 
may not preserve the geometric properties of the system 
accurately [11]. 

Variational integrators provide a systematic method 
of developing geometric numerical integrators for La- 
grangian/Hamiltonian systems [16]. As it is derived from a 
discrete analogue of Hamilton's principle, it preserves sym- 
plecticity and the momentum map, and it exhibits good total 
energy behavior. Lie group methods conserve the structure 
of a Lie group configuration manifold as it updates a group 
element using the group operation [17]. 

These two methods have been unified to obtain a Lie group 
variational integrator for Lagrangian/Hamiltonian systems 
evolving on a Lie group [13]. This preserves symplecticity 
and group structure of those systems concurrently. It has been 
shown that this property is critical for accurate and efficient 
simulations of rigid body dynamics [12]. 

In this section, we develop a Lie group variational inte- 
grator for a 3D elastic string pendulum. We first construct a 
finite element model, and derive an expression for a discrete 
Lagrangian, which is substituted into the discrete-time Euler- 
Lagrange equations on a Lie group. 



S&rb = 



M r{l)- M Rpctl + M Rn^pc - Mges ■ 6r{l)i_ Finite Element Model 



-jn - MpcR^ r{l) + Mpc^R^ r{l) 
-Mpc(lR^r{l) + MgpcR^ez - Cun 



■ V 



V dt, (9) 

where we repeatedly use the property: y-xz = zy- x for any 

x,y,z £ 

From (IT) and (|9]l, the variation of the action integral is 
given by (5© — 5(&str + <5©r&, and it is equal to zero for any 



We discretize the string by N one-dimensional line el- 
ements. Thus, the unstretched length of each element is 
u = j^. A natural coordinate ( G [0, 1] in the a-th element 
is defined by ^ = i(s — u{a — 1)). Let Sq,Si be shape 
functions given by S'o(C) — 1 ^ C> and 5i(C) = C- These 
shape functions are also referred to as tent functions. The 
position vectors for the end nodes of the a-th element are 
given by rk,a, r^^a+i when t — kh for a fixed time step h. 



Using this finite element model, the position vector r(s, t) 
of a material point in the a-th element is approximated as 
follows: 



r{s, t) = So{C)rk,a + Si{C)rkM+i = ffc,a(C)- 



(13) 



Note that rfc,a(0) rk^a and rfc,a(l) = rkM+i- The partial 
derivative with respect to s is given by 

, , , dr(s, t) dC I , . , 

r'i-^^t) = = -{rk,a+i-rk,a)^r',^a- (14) 

The partial derivative with respect to t is approximated by 

r{s,t) = i(5'o(C)Arfc,a + ^i(C)Ar-fc,a+i) =«m(C), (15) 

where the Delta-operator represents a change for a time step, 

i.e. Arfe,a = rfe+i,a - rk^a- 

B. Discrete-Lagrangian 

Using these finite element model, a configuration of the 
discretized 3D elastic pendulum at t ^ kh is described 
by 9k = {vks, ■ ■ ■ ,rk,N+i,Rk), and the coiTesponding 



configuration manifold is G = (] 



S0(3). 



We define a discrete-time kinematics equation as follows. 
Define fk = {Ark,i, . . . , Ark,N+i,Fk) G G for Ark,a G 
and Fk G S0(3) such that gk+i = gkfk and G acts on itself 
by the diagonal action: 

(rk+i,!, ■ ■ ■ , rk+i,N+i, Rk+i) 
= {rk,i + Arfc,i, . . . ,rfc,jv+i + Ark^N+i, RkFk). (16) 

Therefore, fk represents the relative update between two 
integration steps. This ensures that the structure of the Lie 
group configuration manifold is numerically preserved since 
gk is updated by fk using the right Lie group action of G 
on itself. 

A discrete Lagrangian Ld{gk, fk) : G x G ^ R is an 
approximation of the Jacobi solution of the Hamilton-Jacobi 
equation, which is given by the integral of the Lagrangian 
along the exact solution of the Euler-Lagrange equations over 
a single time step: 



Ldigk,fk) 



L{g{t),r'{t)~g{t))dt, 



where g{t) : [0, /i] G satisfies Euler-Lagrange equations 
with boundary conditions g(0) = gk, g{h) = gkfk- The 
resulting discrete-time Lagrangian system, referred to as 
a variational integrator, approximates the Euler-Lagrange 
equations to the same order of accuracy as the discrete 
Lagrangian approximates the Jacobi solution. 

Substituting (fT3]l-(fT5]l into the continuous-time Lagrangian 
given by (|6]l, the contribution of the a-th element to the 
discrete Lagrangian is chosen as follows. 



f^\EA{\\r',^J~V, 



) - fJ-g rk.aiC) ■ ea udC 



This is given by 
1 



mArk.a ■ Ark,a + —mArk,a ■ Ar^ a+i 
6/1 oh 



1 
h 



mArk,a+i ■ ArkM+i 



+ -mg{2rk,a + 'irk^+i + Ark^ + Ark,a+i) ■ 63 

- i/iK(||rfc^a+i - rk^aW - 

- i/i ||rfc,a+i + Ark,a+i ~ rk^ ~ Ark,a\\ - uf ), 



(17) 



where m — /lu, k = So, the contribution of the string 
to the discrete Lagrangian is Ld^ ~ X^^i-^dfca- The 
contribution of the rigid body to the discrete Lagrangian is 
chosen as follows. 

Ld,,^, = ^MArkM+i ■ ArkM+i + ^tr[(/ - Fk)Jd] 
M 

+ —Ark,N+i ■Rk{Fk-I)pc 
h 

+ -i^Mg [rkM+i + RkPc) ■ 63 

+ ^Mg {rk.N+i + ArkM+i + RkFkpc) ■ 63, (18) 

where Jd £ K'^^'' is a nonstandard inertia matrix defined by 
Jd ~ 5tr[J] /3x3 — J, as introduced in [15]. 

From ( [TT] ). ( fTsT i, the discrete Lagrangian of the 3D elastic 
string pendulum is as follows. 

Ldk{9k,fk) = Ld^^t^igk, fk) + Ld^^^igk, fk) 

N 

= X! ''9k, fk) + Ld^ ,^, {gkjk)- (19) 
0=1 

C. Discrete-time Euler-Lagrange Equations 

For a discrete Lagrangian on G x G, the following discrete- 
time Euler-Lagrange equations, referred to as a Lie group 
variational integrator, were developed in [13]. 



Tei-/fc_i ■ D/^._iLdj^_i-Ad*-i ■ (T*L/^ ■'Df^Ld^ 



(20) 



+ T*Lg^ • Dg^Ld^ - 0, 
9k+i = gkfk, (21) 

where TL : TG TG is the tangential map of the left 
translation, D/ represents the derivative with respect to /, 
and Ad* : G x g* g* is co-Ad operator [18]. 

Using this result, we develop a Lie group variational 
integrator for a 3D elastic string pendulum. For / = 
{Ari, . . . ,ArN+i,F) G G and p = {pi, . . . ,pN+i,n) G 
g* ~ (R3)^+^ x K^, the co-Ad operator is given by 
Ad}-ip= {pi,...,pn+i,F7t). 



Derivatives of the discrete Lagrangian: We now obtain 
expressions for the derivatives of the discrete Lagrangian. 
The derivatives of the discrete Lagrangian of the a-th ele- 
ment, given by ( fTTI i. with respect to Ar^ „ and Ar/j q+i are 
given by 



D 



m{Ark,a + ^Arfc,a+i 



-mges 



+ 2 V^fc+l.a, 

= -!-m(Arfc_a+i 



3h 



-Ar 



k,a) 



k+l,a- 



where VV, 



k,a 



\\x\[—u 

k " u u X for X 



Then, from ( fT9] l, the derivative of the discrete Lagrangian 
with respect to Ar^. q, for a € {2, . . . , N}, is given by 

= — m(Arfe.a„i + 4Arfc,a + Ar^.^+i) 
on. 

+ ^mge3 + ^Vy,^i^, - ^Vy^+i^,_i. (23) 

Similarly, the derivative of the discrete Lagrangian with 
respect to rfc ^, for a G {2, . . . , N}, is given by 



(Vy;r,a-l+VVl!Vi,,_i). 



(24) 



Next, we find the derivatives of the discrete Lagrangian 
with respect to Ar^. jv+i and rk^N+i- They are contributed 
by the A^-th string element and the rigid body, and they can 
be obtained from ( fTSI l and (|22] | as follows. 



— i?fc(i^fe-/)p, + -(M + -).ge3 



fc+l,Arj 



(25) 



D 



-VV, 



fc.Af 



(26) 



Now, we find the derivatives of the discrete Lagrangian 
with respect to Fk and Rk- From ( fTSl ), we have 

Di.,irf,. • (5Ffc = itr[-5Ffe Jd] + ^Arfc,jv+i • RkSFkPc 
h 

+ -MgRkSFkPc ■ 63 
= hrl-SFkJd] + Ak ■ SFkPc, 



where ^fe = ^RlArk,N+i + ^MgRle^. The variation of 



Ffe can be written as JFfe = FfcCfe for C,k G K''- Therefore, 
this can be written as 



1 



-FkCkJd 



Ak ■ FkCkPc 



By repeatedly applying the following property of the trace 
operator, tv[AB] = tr[S^] = tr[^^i3^] for any A, B e 
M'^^'^, the first term can be written as tr[~FkCkJd] = 
^hCkJdFk] = triCkF^Jd] = -itr[a(J<iFo - F^J^)]. 
Using the property of the hat map, x'^y ~ — itr[a;y] for any 
x,y e R^, this can be further written as {{JdFk — F'j^ JdY) ■ 



Ck- As y ■ xz 
can be wr 
we obtain 



W^, the second term 

can be written as FjFAk ■ CkPc — PcFk Ak ■ Ck- Using these. 



zy ■ X for any x,y, z £ 



1 



TJt_Ffc • DFfcidfc - J^i-^dFk - F^JdY 



PcF^Ak. 



(27) 



(22) '■^^ CO- Ad operator yields 



Ad;,x 



1 



dk) 



{FkJd - JdFk ) 



FkPcAk- 
(28) 



Similarly, we can derive the derivative of the discrete 
Lagrangian with respect to Rk as follows. 

M 

TJLfl, • T>R,Ld, = —{{Fk - I)perRlArk,N+i 

+ f^Mgp.Rles + ^MgKpcRles. 

(29) 

Discrete-time Eiiler-Lagrange Equations: Substituting 
(l23Tl-(|29]l into (l20]i-(l2ni. we obtain discrete-time Euler- 
Lagrange equations for a 3D elastic string pendulum as 
follows. 



-^m(A^rfc,<j-i 
o/i 



- AA'^rk,a + A'^rk,a+i) 
hVV^.a-i - hVVl, = 0, 



-(M- 



— hmge^ 

m 2 1 A 2 

— )A rk,N+i + -zrmA rk.N 
6 on 



k.N 



-M{RkFk ~ 2Rk + Rk~i)Pc - h{M + -)gez 



-^{FkJd — JdFk — JdFk^l 



Fl.Jdf 



—p.RlAW 



N+i - hMgpcRles 

rk+l.a = rk,a + Arfc,a, 
^fc+l = RkFk- 



(30) 



= 0, 
(31) 

(32) 

(33) 
(34) 



whereAVfe,a = Ar^^a - Art-i.a = rk+i,a-'2rk,a + rk- 



l,a. 



-pu, K = and VVj 



-X for X 



k,a — 

rk,a+i ^ ^k.a- Equation ( [30] l is satisfied for a S {2, . . . , N}, 
and ( |33] | is satisfied for a e {2, . . . , + 1}. For any fc, the 
vector rk.i = since the pivot is fixed. 

For given {gk-i, fk-i)^ 9k is explicitly computed by (|33] l 
and ( l34l i. The update fk is computed by a fixed point 
iteration for Fk'. we select an initial guess of Fk', Ark,a 
is obtained by solving ( |30] | and OTl i. which requires the 
inversion of a fixed 3A^ x 3A^ matrix; a new Fk is computed 
by solving the implicit equation (l32] i: these are repeated until 
Fk converges. When solving the implicit equation ( |32] |. we 
first express Fk on R'^ using the Cayley transform, and apply 
Newton's iteration (See Section 3.3.8 in [13]). These yields 



(a) t e [0, 1.25] 



(b) t G [1.25,2.5] 



(a) Energy transfer (kinetic energy of 
the rigid body: red, kinetic energy of 
the string: black, gravitational poten- 
tial: green, elastic potential: blue) 



(b) Total energy 




(c) t e [2.5,3.75] 



(d) t e [3.75, 5] 



2 3 



Fig. 2. Snapshots of a 3D elastic string pendulum maneuver. Strain energy 
distribution is illustrated by color shading (An animation is available at 
|http://my.fit.edu/'taeyoung) 

a Lagrangian flow map {gk-i,fk-i) ^ {gk,fk), and they 
are repeated. 

V. Numerical Example 

We now demonstrate the computational properties of the 
Lie group variational integrator developed in the previous 
section by considering a numerical example. The material 
properties of the string are chosen to represent a rubber string 
as follows [6]. 

71= 0.025 kg/m, ; = lm, EA = 40N. 

The rigid body is chosen as an elliptic cylinder with a 
semimajor axis 0.06 m, a semiminor axis 0.04 m, and a 
height 0.1m. Its properties are as follows. 

M = 0.1kg, = [0.04, 0.01, 0.05] m, 

0.38 -0.04 -0.20] 
J= -0.04 0.58 -0.05 kgm^ 
-0.20 -0.05 0.30 

Initially, the string is aligned to the horizontal ei axis at rest, 
and the rigid body has an initial velocity [0, 0.2, —0.5] m/s. 
We use iV = 20 elements. Simulation time is T = 5 seconds, 
and time step is h ~ 0.0001 second. 

Fig. [2] illustrates the resulting maneuver of the 3D elastic 
string pendulum. As the point where the string is attached 
to the rigid body is displaced from the center of mass of 
the rigid body, the rigid body dynamics are directly coupled 
to the elastic string dynamics, which yields the illustrated 
complex maneuver. 

Fig. [3] shows the corresponding energy transfer, total 
energy, total angular momentum deviation, orthogonality 



(c) Deviation of the total angular mo- (d) Orthogonality error of rotation 
mentum about the gravity direction matrices ]|/ — R?" R\\ 




(e) Velocity / angular velocity of the (f) Stretched length of the string 
rigid body (second components) 

Fig. 3. Numerical simulation of a 3D elastic string pendulum 



errors of rotation matrices, velocities of the rigid body, and 
the stretched length of the string. As shown in Fig. |3(b)[ the 
computed total energy of the Lie group variational integrator 
oscillates near the initial value, but there is no increasing or 
decreasing drift for long time periods. This is due to the 
fact that the numerical solutions of symplectic numerical 
integrators are exponentially close to the exact solution of 
a perturbed Hamiltonian [19]. The value of the perturbed 
Hamiltonian is preserved in the discrete-time flow. The Lie 
group variational integrator preserves the momentum map as 
in Fig. |3(c)[ and it also preserves the orthogonal structure 
of rotation matrices accurately. The orthogonality errors, 
measured by ||/ — R^R\\, are less than 2 x 10^^'^ in Fig. 

m 

These show that the Lie group variational integrator pre- 
serves the geometric characteristic of the 3D elastic string 
pendulum accurately even for the presented complex ma- 
neuver that has nontrivial energy transfer between different 
dynamic modes. 



VI. Conclusions 

We have developed continuous-time equations of motion 
and a geometric numerical integrator, referred to as a Lie 
group variational integrator, for a 3D elastic string pen- 
dulum. The continuous-time equations of motion provide 
an analytical model that is defined globally on the Lie 
group configuration manifold, and the Lie group variational 
integrator preserves the geometric features of the system, 
thereby yielding a reliable numerical simulation tool for 
complex maneuvers over a long time period. 

These can be extended to include the effects of control 
inputs by using the discrete Lagrange-d'Alembert principle 
[20], which modifies the discrete Hamilton's principle by 
taking into account the virtual work of the external control 
inputs. When applied to an optimal control problem, this 
allows us to find optimal maneuvers accurately and effi- 
ciently, as there is no artificial numerical dissipation induced 
by the computational method. Furthermore, optimal large- 
angle rotational maneuvers can be easily obtained without 
singularities and complexity associated with local parame- 
terizations, since the configuration is represented globally on 
the Lie group [21]. 
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